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INFLUENCE  OF  MICROTURBULENCE  ON  EARLY  TIME  HANE  STRUCTURE 
I.  INTRODUCTION 

An  important  issue  to  the  Defense  Nuclear  Agency  is  understanding  the 
physical  processes  involved  in  the  structuring  of  the  ionosphere  following 
a  HANE.  The  large-scale,  long-lasting  ionization  irregularities  produced 
by  a  HANE  can  adversely  Impact  radar  and  communications  systems  (i.e., 
scintillations).  Thus,  it  is  crucial  to  understand  the  structuring 
mechanisms  associated  with  HANES  in  order  to  obtain  a  predictive  capability 
which  could  aid  the  operation  of  these  systems  in  a  nuclear  environment. 
Recently,  there  has  been  renewed  interest  in  the  DNA  community  concerning 
early  time  HANE  processes  (t  <  few  sec).  For  example,  DNA  is  presently 
funding  a  major  laser  experiment  at  NRL  to  explore  the  physics  associated 
with  early  time  HANE  evolution. 

Among  the  various  early  time  processes  to  be  understood,  the 
occurrence  of  plasma  microinstabilities  and  of  structure  producing 
processes  (i.e.,  interchange  instabilities)  are  two  important  areas.  To 
our  knowledge,  these  two  areas  have  been  treated  somewhat  independently  of 
each  other.  The  studies  of  microturbulence  have  been  directed  at 
explaining  the  coupling  of  debris  and  air  plasmas  (i.e.,  momentum 
exchange),  and  on  the  heating  of  the  plasmas  (i.e.,  energy  exchange)  (Lampe 
et  al.,  1975).  On  the  other  hand,  the  studies  of  early  time  structuring 
processes,  which  have  been  limited  to  date,  do  not  directly  include 
microturbulence  effects.  It  should  be  noted  that  microturbulent  processes 
have  been  included  indirectly  in  an  early  time  structure  analysis.  The 
work  of  Brecht  et  al.  (1982)  addresses  the  growth  of  Rayleigh-Taylor 
instabilities  based  upon  the  conditions  prescribed  by  the  microphysics 
(e.g.,  thickness,  temperature,  density  within  the  coupling  shell). 

The  purpose  of  this  report  is  to  suggest  a  direct  effect  of 
microturbulence  on  the  evolution  of  early  time  Interchange  instabilities. 
Namely,  that  microinstabilities  can  provide  a  dissipation  mechanism  (via 
the  anomalous  transport  properties  associated  with  them)  which  can  cause 
saturation  of  an  interchange  instability  through  a  nonlinear  mode  coupling 
process.  In  order  to  demonstrate  this  point,  we  consider  a  simplified 
model.  We  consider  a  two-fluid  plasma  model  (ions  and  electrons)  which 
Manuscript  approved  January  24,  1984. 


incorporates  a  gravitational  field  (g)  and  a  density  inhomogeneity  (Vn). 
We  develop  a  set  of  nonlinear  mode  coupling  equations  for  this  situation 
and  show  that  the  intermediate  regime  of  the  instability  can  saturate  by 
transferring  energy  from  the  growing  modes  to  the  damped  modes  because  of 
the  nonlinear  electron  E  x  B  motion.  The  crucial  point  in  the  analysis  is 
that  anomalous  electron-ion  collisions  provide  the  dissipation  mechanism 
for  the  damped  modes. 

The  organization  of  the  paper  is  as  follows.  In  the  next  section  we 
derive  the  nonlinear  mode  equations.  In  Section  III  we  present  the  results 
for  a  two  mode  coupling  process.  Finally,  in  Section  IV  we  discuss  the 
implications  of  our  theory  for  early  time  HANE  processes  and  describe  our 
future  efforts  in  this  area. 

II.  DERIVATION  OF  MODE  COUPLING  EQUATIONS 

The  geometry  and  plasma  configuration  used  in  the  analysis  are  shown 
in  Fig.  1.  We  take  the  ambient  magnetic  field  to  be  in  the  z  direction 

A 

(B  *  B  ez),  and  the  density  to  be  inhomogeneous  in  the  x  direction  (n  * 
n(x)).  The  fundamental  assumptions  used  in  the  analysis  are  the 

following.  We  assume  perturbed  quantities  to  be  proportional  to  expfiCl^x 
+  kyy  -  u>t)]  with  k^L  »  1,  kyL  »  1,  k^^  «  1,  kyPj^  «  1,  and  u  «  ^ 

where  L  *  (3£n  n/3x)  *  is  the  density  gradient  scale  length  p.  - 

1/2  X 
(T^/m^)  /fl^  is  the  mean  ion  Larmor  radius,  and  ft^  -  eB/n^c  is  the  ion 

cyclotron  frequency.  We  consider  a  weakly  collisional  plasma  such  that 

v  ,/ft  «  1,  v.  /ft,  «  1,  and  v,,/ft,  «  1  where  v  „  represents  a  collision 
ei  e  ie  i  ii  i  a6 

frequency  between  the  a  and  S  species,  and  ftg  ■  eB/mgc  is  the  electron 
cyclotron  frequency.  Here,  the  collisions  are  attributed  to  wave-particle 
interactions,  i.e.,  anomalous  collisions.  Finally,  we  assume  quasi¬ 
neutrality  (ne  »  n^). 

The  equations  used  in  the  analysis  are  continuity,  momentum  transfer, 
and  charge  conservation: 
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Plasma  configuration  and  slab  geometry.  The  magnetic  field  is 

A 

in  the  z-direction  (B  -  B  e  ] ,  the  density  gradient  and  gravity 
are  in  the  x-direction  ( Vn  -  3n/3x  ex  and  g  -  g  ex)  ,  and  the 
wave  vector  is  in  the  xy  plane  transverse  to  B. 


W  “  Mn(vt  -  Ye)]  -  0 
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1/2 

where  a  denotes  species  (e,i),  ■  (T^/m^)  is  the  thermal  velocity,  and 

is  the  ion  stress  tensor  such  that  (Stix,  1969) 
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The  first  term  in  Eq.  (5)  represents  finite  ion  Larrnor  radius  corrections 
(Roberts  and  Taylor,  1962),  while  the  second  term  represents  ion  viscosity 
effects  (i.e.,  shear  stress).  The  compressional  stress  term  in  Eq.  (5)  has 
been  neglected  since  it  is  smaller  by  a  factor  kp^  than  the  first  two  terms 
(Stix,  1969). 

To  lowest  order  the  equilibrium  particle  drifts  are  given  by  (from 
Eqs.  (2)  and  (3)) 
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de  y 
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where  V.  ■  fv  /8  ]3in  n/3x  is  the  diamagnetic  drift  velocity  of  species  a 
aa  a  a 

and  Vg  -  -  g/8^  is  the  ion  gravitational  drift  velocity. 

We  now  perturb  Eqs.  (1)  -  (4)  about  this  equilibrium  and  let  n  -  nQ  + 
Sn,  ,  and  V  •  V  +  6V  .  +  6V  «.  Here,  SV  ,  is  the  first  order 

velocity  perturbation  and  dV^  *s  t*ie  secon{*  order  velocity  perturbation 
(i.e.,  proportional  to  small  factors  such  as  ve^/ne»  vie^i»  vii^i»  ^pi* 

8^  ^d/dt).  The  first  order  perturbed  velocities  are 
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where  ip  ■  edij>/T,  +  5n/nn. 


The  second  order  velocities  are 
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Substituting  Eqs.  (8)  -  (11)  into  the  electron  continuity  (Eq.  (1))  and 
charge  conservation  equations  (Eq.  (4)),  we  arrive  at  the  coupled  set  of 
equations  for  6n  and  6$ 
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where  D_  ■  v  ,p”  is  the  electron  diffusion  coefficient  and  p  ■  fv,/ 
e  el  es  es  v  i 

Vne)(l  +  Te/Ti)* 

Thus,  Eqs.  (12)  and  (13)  represent  nonlinear  mode  coupling  equations 
for  the  Rayleigh-Taylor  instability  in  the  intermediate  wavelength  regime 
(i.e.,  kL  »  1  and  kpj^  «  1).  The  physics  contained  in  Eqs.  (12)  and  (13) 
is  described  as  follows.  The  second  term  on  the  LHS  of  Eq.  (12)  leads  to 
growth  of  the  density  perturbation  while  the  third  term  leads  to  damping 
because  of  electron  diffusion.  The  term  on  the  RHS  of  Eq.  (12)  is  the 
electron  E  x  B  nonlinearity.  The  first  term  on  the  LHS  of  Eq.  (13) 
represents  the  differential  motion  of  the  electrons  and  ions  in  the  y 
direction  because  of  gravity.  The  first  term  in  brackets  is  the  ion 
polarization  term  while  the  second  term  contains  ion  viscosity  effects. 


III.  ANALYSIS  OF  MODE  COUPLING  EQUATIONS 
A.  Linear  Theory 

We  linearize  Eqs.  (12)  and  (13)  by  neglecting  the  electron  E  x  B 
nonlinearity  in  Eq.  (12).  We  assume  that  perturbations  are  proportional  to 

exp [ i (Ic^x  +  kyy  -  tot)].  From  Eq.  (12),  it  can  be  shown  that 

0  to  +  iv.kp  i 
ei  es 

where  L  -  (Vn^/n^)  ^  »  (3  in  ng/3x)  *  is  the  scale  length  of  the  density 

gradient.  Similarly,  from  Eq.  (13)  it  can  be  shown  that 

-  [»g  +  +  «g  +  i  0-3  v^pj)]  p-  «  (15) 

k2p^  [o>  +  uig  +  i  0.3  viife2P^] 

where  u  *  k  g/fl..  Combining  Eqs.  (14)  and  (15),  we  obtain  the  linear 
g  y  i 

dispersion  equation 

u>2  +  fau  +  c  -  0  (16) 


where 


+  °-3 

(17) 

*  t2  2i 

L  “  iveik  pesJ  x 

(18) 

.2  2-i 

iik  PJ 

$  t  ‘ [  V 


and  vi0  -  Vdi  +  vg,  Vdl  -  /fJ^L  and  Vg  ■  -  g/fl^.  Equation  (16)  has  the 
solution 


-  b/2  +  (b2/4  -  c) 


which  is  rather  messy  given  the  complexity  of  b  and  c  (Eqs.  (17)  and 
(18)).  In  order  to  gain  insight  into  the  solution  of  Eq.  (19),  it  is  best 
to  consider  several  limiting  cases. 


1.  Collisionless  plasma  (vg^  =*  *  0) 

In  the  collisionless  limit,  the  dispersion  equation  is  given  by 


2 

oo 


-  k  V,n  hi 
y  iO 
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(20) 


which  is  a  well  known  result  (Roberts  and  Taylor,  1962). 
Eq.  (20)  is 


1  i  ,  ,  K  i/2 
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In  the  limit  k  «  4g/L,  the  eigenf requency  is 


CO 


y  r  S'! 1/2 


The  solution  to 


(21) 


(22) 


so  that  instability  results  if  g  and  Vn  are  oppositely  directed,  i.e.,  g/L 

2  2 

<  0.  On  the  other  hand,  in  the  limit  k  >  4g/L,  the  mode  becomes  stable 
(y  ■  0)  and  propagating  (uir  ~  kyVi0)*  Since  it  is  generally  the  case  that 
V a  -  (p^/L)v^,  the  stabilization  criterion  is  roughly 

k2pj  >  4gL/v2  (23) 


and  is  therefore  known  as  'finite  Larmor  radius  stabilization'.  An 
important  point  is  that  the  stable  modes  have  y  m  0  and  are  not  damped, 
i.e.,  Y  <  0.  This  is  a  crucial  consideration  in  the  mode  coupling 
saturation  mechanism  since  damped  modes  are  required  as  energy  sinks. 


2.  Collisional  plasma  (v@^  *  0  and  *  0) 

2  2 

We  now  consider  the  strong  collision  limit  given  by  v^k  p.^ 
kyV10.  For  this  case  the  dispersion  equation 


1»  0.3  »uk2pJ 


-l  f  +  0.3 

*_L  Li 


.2  2  .22 
v  .  .k  p ,v  .k  p 
ii  i  ei  es 


»  10  , 


(24) 


7 


which  has  the  solution 
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Again  instability  can  result  for  g/L  <  0  (oppositely  directed  g  and  Vn)  as 
long  as 


y  I g/L I 

2  2  “> 
k  0.3  v^kpj 


In  the  opposite  limit, 
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the  modes  are  damped,  i.e.,  y  <  0.  The  important  point  is  that  the 

electron-ion  collision  term  provides  a  dissipation  mechanism  and  can 

produce  modes  with  a  negative  growth  rate.  Hence,  anomalous  electron-ion 

collisions  can  produce  damped  modes  in  the  collisional  limit  (y  <  0).  This 

is  in  contrast  to  collisionless  finite  Larmor  radius  stabilization  which 

produces  stable  propagating  waves  (y  *  0).  Finally,  we  note  the  similarity 

of  Eq.  (25)  to  that  of  the  usual  collisional  Rayleigh-Taylor  mode  which  has 

2  2 

a  growth  rate  y  *  -  (k  /k  )(g/Lv.  )•  In  the  present  case,  the 

2  2  y  in 

quantity  0.3  v^k  p^  replaces  (Hudson  and  Kennel,  1975). 

B .  Nonlinear  Theory 

To  illustrate  the  nonlinear  saturation  of  the  Rayleigh-Taylor 
instability  we  consider  the  collisional  limit  and  follow  the  analysis 
prescribed  in  Rognlein  and  Weinstock  (1974).  Mathematically,  this  analysis 
closely  follows  that  of  Chaturvedi  and  Ossakow  (1977),  but  the  underlying 
physical  concepts  and  applications  are  substantially  different.  The 
equations  of  interest  are 
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where  Eq.  (28)  is  the  same  as  Eq.  (12)  and  Eq.  (29)  follows  from  Eq.  (13) 
by  neglecting  ion  inertial  effects. 


Equation  (28)  is  linear  and  one  can  obtain  a  simple  relationship 
between  <Sq>  and  <5n  based  upon  it,  namely, 
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where 
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and  we  have  assumed  kyL  > 

Equation  (27)  is  nonlinear  and  can  be  rewritten  as 
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where  y  is  the  net  linear  growth  (or  damping)  of  modes.  The  term  y  results 
from  a  combination  of  the  second  term  on  the  LHS  of  Eq.  (28)  (the  growing 
term)  and  the  third  term  on  the  LHS  of  Eq.  (28)  (the  damping  term),  and  is 
given  by 

2 

Ig/Ll  _  ,.2.2 


k 


T?  '  veik  pes 


(33) 


k  0,3  Viik  Pi 

for  g/L  <  0. 

We  substitute  a  two-dimensional  perturbation  of  the  form 


—  ■  Altl  sin(kyy  _  »t)  cos  kxx 


(34) 


and  (from  Eq.  (29)) 


1^-  ■  T  A.  .  cos  (k  y  -  oat)  cos  k  x 
A  ^  * » *■  y  x 


(35) 


into  the  nonlinear  term  of  Eq.  (32).  This  procedure  yields  the  following 
result 


CvCvw.'- 


.  9  »  • 
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(36) 


~  TTT  .2  kxAi,i 
0,6  viik  pi  k 


sin  2  k  x, 
x 


that  is,  a  second  spatial  harmonic  in  the  x-direction  (along  the  density 
gradient)  is  produced.  Thus,  we  consider  a  general  perturbation  of  the 
form 


—  -  A.  sin  (k  y  -  ait)  cos  k  x  +  A,  n  sin  2  k  x  (37) 
tiQ  i»i  y  x  ^ ”  x 

Substituting  Eq.  (37)  into  Eq.  (32)  and  making  use  of  Eq.  (30),  we 
arrive  at  the  following  set  of  coupled  nonlinear  differential  equations  for 
the  mode  amplitudes 


at  Ai,i  "  Yi,iAi,i  "  a  ai,ia2,o 


(38) 


and 


—  A 

at  2,o 


T2,0A2,0 


2 
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(39) 


where  the  coupling  coefficient  is 


a 


_y _ 

k2  0.3 


_S_ 


.2  2 

Viik  Pi 


(40) 


The  mode  A^  q  is  linearly  damped  (i.e.,  Y2  q  <  0)  which  is  clear  from 
Eq.  (33)  since  ky  ■  0  for  this  mode.  Thus,  the  action  of  the  nonlinear 
E  x  B  term  is  to  transfer  energy  from  a  growing  mode  (A.  1 )  to  a  damped 

*  f  * 

mode  (A2  q)  via  a  mode  coupling  interaction.  This  is  shown  schematically 
in  Fig.  2.  It  is  possible  then  to  set  up  a  steady  state  situation  such 
that  3A^  ^/3t  «  0  and  3A^  q/3  t  -  0.  One  can  show  that  this  results  in  the 
following  set  of  steady  state  amplitudes 


(42) 
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We  point  out  that  for  Iy^  >  Iyj  qI»  one  finds  that  A2  0  >  Al  i  30 
that  the  damped  modes  can  be  nonlinearily  driven  to  larger  amplitudes  than 
the  linear  driven  modes.  These  results  are  similar  to  those  of  Rognlein 
and  Weinstock  (1974),  Chaturvedi  and  Ossakow  (1977)  and  Chaturvedi  and 
Ossakow  (1979).  Finally,  although  a  rather  simple  two-mode  analysis  has 
been  presented,  numerical  simulations  using  many  modes  for  both  the 
collisional  Rayleigh-Taylor  and  E  x  B  instability  are  consistent  with  this 
analysis  (Keskinen  et  al.,  1980a, b). 


ky 

Figure  2:  Schematic  of  mode  coupling  process.  The  growing  and  damped 
modes  are  shown  in  k  space.  The  coupling  of  the  growing  (1,1) 
mode  and  the  damped  (0,2)  mode  is  also  shown. 
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IV.  DISCUSSION 


The  purpose  of  this  paper  has  been  to  suggest  an  important  influence 
of  microturbulence  on  the  early_  time  structuring  of  a  HANE.  We  consider 
the  question,  how  do  anomalous  collisions  produced  by  microturbulence 
affect  the  development  of  interchange  instabilities?  To  answer  this 
question,  we  have  adopted  a  simple  two-fluid  model  of  the  intermediate 
wavelength  Rayleigh-Taylor  instability  as  a  test  problem.  Specifically,  we 
have  developed  a  two-dimensional  nonlinear  mode  coupling  theory  of  the 
Rayleigh-Taylor  instability  in  the  presence  of  microturbulence.  The 
primary  effect  of  the  microturbulence  is  to  generate  anomalous  collisions 
between  electrons  and  ions,  as  well  as  ions  and  ions;  an  effect  considered 
to  be  very  important  in  the  early  time  evolution  of  a  HANE.  Examples  of 
Instabilities  investigated  for  early  time  HANE  applications  are  the  ion-ion 
counter-streaming  Instabilities,  the  modif ied-two-stream  instability,  the 
ion  acoustic  instability,  and  the  beam  cyclotron  instability  (Lampe  et  al., 
1975).  Lampe  et  al.  (1975)  have  shown  that  these  instabilities  can  produce 
significant  anomalous  collisions  between  particles,  e.g.,  vgi  <  w for  the 

modified  two  stream  and  for  the  magnetized  ion-ion  counter¬ 

streaming  instabilities  where  Upg  is  the  electron  plasma  frequency  and 
is  the  lower  hybrid  frequency.  We  note  that  in  the  model  considered,  the 
concept  of  anomalous  ion-ion  collisions  for  a  single  ion  fluid  is  unclear, 
l.e.,  anomalous  collisionless  viscosity  is  not  well  understood,  but  it  is  a 
valid  concept  for  a  multi-ion  species  plasma.  Thus,  as  noted  in  the 

introduction,  the  problem  addressed  in  this  paper  should  be  considered  as 
an  idealized  test  problem  which  highlights  a  potentially  important  effect 
in  early  time  HANE  physics. 

The  major  result  of  this  work  is  that  the  intermediate  wavelength 
modes  of  the  Rayleigh-Taylor  Instability  can  be  stabilized  via  a  mode 
coupling  process  in  the  collisional  regime.  Basically,  wave  energy  is 
transferred  nonlinear ly  from  growing  modes  to  damped  modes,  and  a  steady 
state  can  result.  Anomalous  electron-ion  collisions  play  a  crucial  role  in 
that  they  provide  a  dissipation  mechanism  which  produces  a  set  of  damped 
modes.  This  is  in  contrast  to  stabilization  of  the  instability  due  to 
finite  Larmor  radius  effects  which  does  not  produce  damped  modes,  but 

rather  propagating  modes  with  y  ■  0. 


A  considerable  amount  of  work  remains  to  be  done  on  this  problem. 
Within  the  context  of  the  simple  model  developed  in  this  paper  two  avenues 
of  research  are  presently  being  pursued.  First,  numerical  analysis  is 
underway  to  consider  the  many-mode  situation  which  is  more  realistic.  In 
this  analysis,  a  pseudo-spectral  code  is  being  used  which  has  proved 
valuable  in  understanding  the  mode  coupling  stabilization  of  the  lower- 
hybrid-drift  instability  (Drake  et  al. ,  1983a, b).  Second,  the  effects  of 
ion  inertia  need  to  be  included  self-consistently  to  understand  the 
transition  from  the  collisionless  to  collisional  regime.  Again,  the 
numerical  techniques  being  used  allow  for  this. 

While  the  above  described  work  will  provide  a  better  understanding  of 
the  nonlinear  stabilization  of  the  Rayleigh-Taylor  instability,  it  does  not 
directly  apply  to  early  time  HANE  phenomena.  It  does  provide  insight  into 
potentially  important  effects,  and  the  development  of  relevant  numerical 
tools.  However,  considerably  different  models  need  to  be  developed  for 
early  time  HANE  behavior.  The  types  of  issues  that  should  be  addressed  are 
as  follows: 

1.  A  stability  analysis  of  early  time  HANE  has  been  performed  by 
Brecht  and  Papadopoulos  (1979).  They  considered  the  development  of  a 
Rayleigh-Taylor  instability  driven  by  forces  associated  with  laminar-like 
acceleration  of  ions  within  the  coupling  shell  and  with  the  curvature  of 
the  magnetic  field  lines.  Another  effect  not  considered  in  Brecht  and 
Papadopoulos  (1979)  is  the  deceleration  of  the  coupling  shell  due  to  line 
tying  currents  to  the  conducting  ionosphere  (Pilipp,  1971;  Fedder,  1980) 
which  could  also  lead  to  a  Rayleigh-Taylor  instability.  The  point  is  that 
there  are  several  mechanisms  that  could  produce  early  time  structure 
instabilities  and  that  a  nonlinear  analysis  of  these  instabilities, 
incorporating  microturbulence  via  anomalous  collisions,  is  warranted. 
However,  this  is  somewhat  more  complicated  than  the  present  calculation 
since  at  least  two  ion  species  are  required  in  the  analysis,  and  electro¬ 
magnetic  effects  may  need  to  be  included. 

2.  Another  possibility  for  the  generation  of  structure  (i.e.,  density 
irregularities)  in  early  time  HANE  is  the  onset  of  the  Richtmeyer-Meshkov 
instability  (Richtmeyer,  1960;  Andronov  et  al. ,  1976).  This  instability 
grows  linearly  with  time,  as  opposed  to  exponentially,  and  can  occur  when  a 
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shock  wave  passes  through  an  Irregular  (corrugated)  interface  from  a  less 

dense  to  a  more  dense  fluid.  The  turbulence  generated  from  this 

instability  may  also  be  Influenced  by  microturbulence,  as  in  the  case  of 

■  ^ 

the  more  standard  interchange  instability  case. 

3.  It  has  also  been  suggested  that  the  weapon  debris  is  'pre¬ 
structured'  by  the  'nuts  and  bolts'  of  the  device  itself  (Chesnut,  private 
communication).  The  subsequent  evolution  of  this  structure  could  also  be 
affected  by  microturbulence.  For  example,  the  wavenumber  spectrum  of  the 
turbulence  could  be  altered  by  the  microturbulence  through  a  mode  coupling 
process. 

In  conclusion,  we  have  shown  that  microturbulence,  via  its  associated 
anomalous  collisional  effects,  can  influence  the  nonlinear  development  of 
the  collisional  Rayleigh-Taylor  instability.  Namely,  it  can  lead  to  the 
stabilization  of  the  intermediate  wavelength  modes  through  a  two- 
dimensional  mode  coupling  process.  Anomalous  electron-ion  collisions 
provide  the  needed  dissipation  mechanism.  Thus,  microturbulence  could  play 
an  important  role  in  the  development  of  structure  for  early  time  HANE 
phenomena.  We  have  given  some  examples  of  the  structure  problems  that 
could  be  affected  by  microturbulence  and  are  relevant  to  early  time  HANE 
processes.  Admittedly,  some  of  these  ideas  are  speculative  at  this  time, 
but  it  is  clear  that  further  research  is  needed  in  this  area. 
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AEROSPACE  DEFENSE  COMMAND/XPD 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AFB,  CO  80912 
OICY  ATTN  XPDQQ 
OICY  ATTN  XP 


a 


AIR  FORCE  GEOPHYSICS  LABORATORY 


HANSCOM  AFB,  MA  01731 


OICY  ATTN  OPR  HAROLD  GARDNER 
OICY  ATTN  LKB  KENNETH  S.W.  CHAMPION 
OICY  ATTN  OPR  ALVA  T.  STAIR 
OICY  ATTN  PHD  JURGEN  BUCHAU 
OICY  ATTN  PHD  JOHN  P.  MULLEN 


AF  WEAPONS  LABORATORY 
KIRTLAND  AFT,  tiM  87117 
OICY  ATTN  SUL 

OICY  ATTN  CA  ARTHUR  H.  GUENTHER 
OICY  ATTN  HTYCE  1LT.  G.  KRAJEI 


AFTAC 

PATRICK  AFB,  FL  32925 
OICY  ATTN  TF/MAJ  WILEY 
OICY  ATTN  TN 


AIR  FORCE  AVIONICS  LABORATORY 
WRICHT-PATTERSON  AFB,  OH  45433 
OICY  ATTN  AAD  WADE  HUNT 
OICY  ATTN  AAD  ALLEN  JOHNSON 


DEPUTY  CHIEF  OF  STAFF 
RESEARCH,  DEVELOPMENT,  &  ACQ 
DEPARTMENT  OF  THE  AIR  FORCE 
WASHINGTON,  D.C.  20330 
OICY  ATTN  AFRDQ 


HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731 
OICY  ATTN  J.  DEAS 


HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION/YSEA 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01732 
OICY  ATTN  YSEA 


V 


J»! 


HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION/DC 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731 

01CY  ATTN  DC  ICC  MAJ  J.C.  CLARK 

( 

COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRICHT-PATTERSON  APB,  OH  45433 
OICY  ATTN  NICD  LIBRARY 

01CY  ATTN  ETDP  B.  BALLARD 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
CRIFFISS  AFB,  NY  13441 

OICY  ATTN  DOC  LIBRARY/TSLD 
OICY  ATTN  OCSE  V.  COYNE 

SAMSO/SZ 

POST  OFFICE  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANCELES,  CA  90009 
(SPACE  DEFENSE  SYSTEMS) 

OICY  ATTN  SZJ 

STRATEGIC  AIR  COMMAND/XPFS 
OFFUTT  AFB,  NB  68113 

OICY  ATTN  ADWATE  MAJ  BRUCE  BAUER 
OICY  ATTN  HRT 

OICY  ATTN  DOK  CHIEF  SCIENTIST 

SAMSO/SK 
P.O.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

OICY  ATTN  SKA  (SPACE  COMM  SYSTEMS) 
M.  CLAVIN 

SAMSO/MN 

NORTON  AFB,  CA  92409 
(MINUTEMAN) 

OICY  ATTN  MNNL 

COMMANDER 

XOME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 

OICY  ATTN  EBP  A.  LORENTZEN 

DEPARTMENT  OF  ENERCY 
LIBRARY  ROOM  G-042 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  FOR  A.  LABOWITZ 


DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.O.  BOX  5400 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR  D.  SHERWOOD 

EG&G,  INC. 

LOS  ALAMOS  DIVISION 

P.O.  BOX  809 

LOS  ALAMOS,  NM  85544 

OICY  ATTN  DOC  CON  FOR  J.  BREEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  TECH  INFO  DEPT 
OICY  ATTN  DOC  CON  FOR  L-389  R.  OTT 
OICY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 
OICY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS,  NM  87545 

OICY  ATTN  DOC  CON  FOR  J.  WOLCOTT 
OICY  ATTN  DOC  CON  FOR  R.F.  TASCHEK 
OICY  ATTN  DOC  CON  FOR  E.  JONES 
OICY  ATTN  DOC  CON  FOR  J.  MALIK 
OICY  ATTN  DOC  CON  FOR  R.  JEFFRIES 
OICY  ATTN  DOC  CON  FOR  J.  ZINN 
OICY  ATTN  DOC  CON  FOR  ?.  KEATON 
OICY  ATTN  DOC  CON  FOR  D.  WESTERVELT 
OICY  ATTN  D.  SAPPEHFIELD 

SANDIA  LABORATORIES 
P.O.  BOX  5800 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR  W.  BROWN 

OICY  ATTN  DOC  CON  FOR  A.  THORNBROUCH 

OICY  ATTN  DOC  CON  FOR  T.  WRIGHT 

OICY  ATTN  DOC  CON  FOR  D.  DAHLGREN 

OICY  ATTN  DOC  CON  FOR  3141 

OICY  ATTN  DOC  CON  FOR  SPACE  PROJECT  DIV 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  B.  MURPHEY 
OICY  ATTN  DOC  CON  FOR  T.  COOK 


OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  DR.  TO  SONG 
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OTHEit  GOVERNMENT 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  D.C.  20234 

OICY  (ALL  CORRES:  ATTN  SEC  OFFICER  FOR) 

INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  &  INFO  ADMIN 
BOULDER,  CO  S0303 

01CY  ATTN  A.  JEAN  (UNCLASS  ONLY) 

01CY  ATTN  W.  UTLAUT 
01CY  ATTN  D.  CROHBIE 
Old  ATTN  L.  BERRY 

NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 
01CY  ATTN  R.  GRUBB 
Old  ATTN  AERONOMY  LAB  G.  REID 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
Old  ATTN  I.  CARFUNKEL 
01CY  ATTN  T.  SALMI 
Old  ATTN  V.  JOSEPHSON 
Old  ATTN  S.  BOWER 
Old  ATTN  D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

Old  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 
Old  ATTN  L.  SLOAN 
01CY  ATTN  R.  THOMPSON 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  94701 
OICY  ATTN  J.  WORKMAN 
01CY  ATTN  C.  PRETTIE 
OICY  ATTN  S.  BRECHT 


BOEING  COMPANY,  THE 
P.O.  BOX  3707 
SEATTLE,  WA  98124 
OICY  ATTN  G.  KEISTER 
OICY  ATTN  D.  MURRAY 
OICY  ATTN  G.  HALL 
OICY  ATTN  J.  KENNEY 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
OICY  ATTN  D.3.  COX 
OICY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20734 
Old  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

OICY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
Old  ATTN  H.  LOCSTON 
OICY  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  Wilshlre  Blvd. 

Santa  Monica,  Calif  90401 
OICY  ATTN  C.8.  GABBARD 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 

OICY  ATTN  M.H.  30RTNER  SPACE  SCI  LAB 

CENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 
OICY  ATTN  F.  REIBERT 
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GENERAL  ELECTRIC  TECH  SERVICES  CO.,  INC. 
HMES 

COURT  STREET 
SYRACUSE,  NY  13201 
01CY  ATTN  G.  MILLMAN 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
01CY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
01CY  ATTN  TECHNICAL  LIBRARY 
OICY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYL VANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

OICY  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

OICY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

OICY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  22311 
OICY  ATTN  J.M.  AEIN 
OICY  ATTN  ERNEST  BAUER 
OICY  ATTN  HANS  WOLFARD 
OICY  ATTN  JOEL  BENGSTON 

INTL  TEL  &  TELECRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY.  NJ  07110 

OICY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TORREYANA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO,  CA  92138 

OICY  ATTN  J.L.  SPERLING 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

OICY  ATTN  DOCUMENT  LIBRARIAN 
OICY  ATTN  THOMAS  POTEMRA 
OICY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINCS,  CO  80933 
OICY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
OICY  ATTN  DAS I AC 
OICY  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMARA 
OICY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
OICY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-12 
OICY  ATTN  W.L.  IMHOF  DEPT  52-12 
OICY  ATTN  RICHARD  C.  JOHNSON  DEPT  52-1 
OICY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 
OICY  ATTN  R.  HEFFNER 


M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

OICY  ATTN  DAVID  M.  TOWLE 
OICY  ATTN  L.  LOUCHLIS 
OICY  ATTN  D.  CLARK 


MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
OICY  ATTN  N.  HARRIS 
OICY  ATTN  J.  MOULE 
01CT  ATTN  GEORGE  MROZ 
OICY  ATTN  W.  OLSON 
01CY  ATTN  R.W.  HALPRIN 
01CY  ATTN  TECHNICAL  LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
01CY  ATTN  P.  FISCHER 
01CY  ATTN  W.F.  CREVIER 
01CY  ATTN  STEVEN  L.  GUTSCHE 
01CY  ATTN  R.  BOGUSCH 
01CY  ATTN  R.  HENDRICK 
01CY  ATTN  RALPH  KILB 
01CY  ATTN  DAVE  SOWLE 
01CY  ATTN  F.  FAJEN 
01CY  ATTN  M.  SCHEIBE 
01CY  ATTN  CONRAD  L.  LONCMIRE 
01CY  ATTN  B.  WHITE 

MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 

ALBUQUERQUE,  NEW  MEXICO  87106 
OICY  R.  STELLINCWERF 
01CY  M.  ALME 
OICY  L.  WRIGHT 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

OICY  ATTN  JOHN  MORGANSTERN 
OICY  ATTN  G.  HARDING 
OICY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
OICY  ATTN  W.  HALL 
OICY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

OICY  ATTN  E.C.  FIELD,  JR. 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMSTRICS ,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01301 

OICY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

OICY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  &  D  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
OICY  ATTN  FORREST  GILMORE 
OICY  ATTN  WILLIAM  B.  WRIGHT,  JR. 
OICY  ATTN  ROBERT  F.  LELEVIER 
OICY  ATTN  WILLIAM  J.  KARZAS 
OICY  ATTN  H.  ORY 
OICY  ATTN  C.  MACDONALD 
OICY  ATTN  R.  TURCO 
OICY  ATTN  L.  DeRAND 
OICY  ATTN  W.  TSAI 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
OICY  ATTN  CULLEN  CRAIN 
OICY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

OICY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

OICY  ATTN  VINCE  TRAPANI 


SCIENCE  APPLICATIONS,  INC. 

1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

OICY  ATTN  LEWIS  M.  LINSON 
01CY  ATTN  DANIEL  A.  HAMLIN 
OICY  ATTN  E.  FRIEMAN 
OICY  ATTN  E.A.  STRAKER 
OICY  ATTN  CURTIS  A.  SMITH 
01CY  ATTN  JACK  MCDOUGALL 

SCIENCE  APPLICATIONS,  INC 
1710  COODRIDGE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 

OICY  ATTN  DONALD  NEILSON 
OICY  ATTN  ALAN  BURNS 
OICY  ATTN  G.  SMITH 
OICY  ATTN  R.  TSUNODA 
OICY  ATTN  DAVID  A.  JOHNSON 
OICY  ATTN  WALTER  G.  CHESNUT 
OICY  ATTN  CHARLES  L.  RINO 
OICY  ATTN  WALTER  JAYE 
OICY  ATTN  J.  VICKREY 
OICY  ATTN  RAY  L.  LEADABRANO 
OICY  ATTN  G.  CARPENTER 
OICY  ATTN  G.  PRICE 
OICY  ATTN  R.  LIVINGSTON 
OICY  ATTN  V.  GONZALES 
OICY  ATTN  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

OICY  ATTN  W.P.  BOQUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA,  CA  93111 
OICY  ATTN  JOHN  ISE,  JR. 

OICY  ATTN  JOEL  GARBARINO 


TRW  DEFENSE  &  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
OICY  ATTN  R.  K.  PLEBUCH 
OICY  ATTN  S.  ALTSCHULER 
OICY  ATTN  D.  DEE 
OICY  ATTN  D /  STOCKWELL 
SNTF/1575 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
OICY  ATTN  W.  REIDY 
OICY  ATTN  J.  CARPENTER 
OICY  ATTN  C.  HUMPHREY 


